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The motion of a pair of counter-rotating point vortices placed in a uniform flow 
around a circular cylinder forms a rich nonlinear system that is often used to model 
vortex shedding. The phase portrait of the Hamiltonian governing the dynamics of 
a vortex pair that moves symmetrically with respect to the centerline — a case that 
can be realized experimentally by placing a splitter plate in the center plane — is 
presented. The analysis provides new insights and reveals novel dynamical features 
of the system, such as a nilpotent saddle point at infinity whose homoclinic orbits 
define the region of nonlinear stability of the so-called Foppl equilibrium. It is pointed 
out that a vortex pair properly placed downstream can overcome the cylinder and 
move off to infinity upstream. In addition, the nonlinear dynamics resulting from 
antisymmetric perturbations of the Foppl equilibrium is studied and its relevance to 
vortex shedding discussed. 
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I. INTRODUCTION 



Flow around a circular cylinder is a classical topic in hydrodynamics that is of funda- 
mental importance to many scientific fields with numerous applications^. Of particular 
interest is the formation, at moderate Reynolds numbers, of vortex eddies behind a circular 
cylinder, which then go unstable at higher Reynolds numbers and evolve into a Karman 
vortex street^. Since an analytic treatment of the problem in terms of the Navier-Stokes 
equation is difficult and the computational cost of direct numerical simulation very high, a 
particularly useful approach to study the basic features of vortex shedding from bluff bodies 
is to consider the dynamics of point vortices in an inviscid fluid. 

A point- vortex model for the formation of two recirculating, symmetric eddies in the wake 
of a circular cylinder was first introduced by Foppl^. He obtained stationary solutions for a 
pair of vortices behind the cylinder in a uniform stream and found that the centers of the 
eddies observed in the experiments lie on the locus of such equilibria — now called the Foppl 
curve. In addition, Foppl found that these equlibria, although stable against perturbations 
that are symmetric with respect to the centerline, were unstable against nonsymmetric per- 
turbations. This instability is believed to constitute the origin of the vortex shedding process 
that leads to the formation of the Karman vortex street^. It was later found out indepen- 
dently by several authors^- that Foppl's stability analysis for symmetric perturbations was 
in error in that the stationary solution behind the cylinder is not exponentially but only 
marginally stable. Physically, marginal stability implies, for instance, that if a splitter plate 
is placed behind the cylinder in the center plane of the wake to suppress vortex shedding^— , 
oscillating forces on the cylinder may still arise owing to the cyclic motion of the vortices 
around their equilibrium position^. 

Despite many contributions to the problem, it is fair to say that the nonlinear dynamics 
of the Foppl system is not yet fully understood. In particular, a more complete picture of 
vortex-pair dynamics in the presence of symmetric perturbations is lacking, and several as- 
pects of the nonlinear dynamics for nonsymmetric perturbations remain unclear. To address 
these two issues is the main motivation of the present paper. It should be emphasized at the 
outset that a better understanding of the dynamical structure underlying the Foppl model 
is of interest not only because of its practical relevance for vortex shedding, but also in its 
own theoretical right from the viewpoint of nonlinear dynamics. 

2 



The Foppl model has inspired a number of studies on several related problems, such 
as the modeling of vortex wake behind slender bodies in terms of multiple pairs of point 
vortices^ - — , the Hamiltonian structure of a circular cylinder interacting dynamically with 
point vortices^— , the control of vortex shedding^ - — , and the stability of symmetric and 
asymmetric vortex pairs over three-dimensional slender conical bodies^ 1 ^ 1 ^. The related 
problem of desingularization of the Foppl pair in terms of vortex patches of finite area 
was also studied^^l A recent review on vortex motion past solid bodies with additional 
references to the Foppl model and related problems can be found in Ref. [28]. 

After formulating the problem of a pair of counter-rotating point vortices placed in a uni- 
form stream around a circular cylinder in Sec. [TTJ we begin our analysis of the Foppl system in 
Sec. IHIl by studying its Hamiltonian dynamics restricted to the invariant subspace where the 
vortices move symmetrically with respect to the centerline. A phase portrait of the system is 
presented that fully characterizes the dynamics within this symmetric subspace. In particu- 
lar, we point out that in addition to the two previously known sets of equilibira, namely, the 
Foppl equilibrium and the equilibrium on the axis bisecting the cylinder perpendicularly to 
the uniform flow, the system possesses a hitherto unnoticed nilpotent saddle at infinity. We 
show furthermore that the homoclinic orbits associated with this nilpotent saddle delimit 
the region of closed orbits around the Foppl equilibrium. We proceed in Sec. IIVI to study 
the linear and nonlinear dynamics resulting from antisymmetric perturbations of the Foppl 
equilibrium. In the linear regime, a mistake that went undetected in Foppl's expressions^ 
for the corresponding eigenvalues is now corrected. As for the nonlinear dynamics, the un- 
stable manifold associated with the Foppl equilibrium is computed numerically and its close 
relation to the vortex shedding instability is pointed out. The linear stability analysis of the 
equilibria on the normal line with respect to symmetric and antisymmetric perturbations 
is also presented — for the first time, it seems — and the respective nonlinear dynamics is in- 
vestigated numerically. A discussion of the physical relevance of our findings and our main 
conclusions are presented in Sec. [V] 

II. PROBLEM FORMULATION 

We consider the motion of a pair of point vortices of same strength and opposite polarities 
around a circular cylinder of radius a and in the presence of a uniform stream of velocity U, 
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FIG. 1. A pair of vortices behind a circular cylinder in a uniform stream. 



as illustrated in Fig. [TJ It is convenient to work in the complex z-plane, where z = x + iy, 
and place the center of the cylinder at the origin. The upper and lower vortices are located 
at positions Z\ = x\ + iyi and z 2 = x 2 + iy 2 , respectively. The complex potential w(z) = 
4>(x, y) + iip(x,y), with <fi being the velocity potential and ip the stream function, is given 
by22 

where T is the circulation of the vortex at Z\ and bar denotes complex conjugation. In Eq. (0Q), 
the first two terms represent the incoming flow and its image (a doublet at the origin) with 
respect to the cylinder, the third term gives the contributions to the complex potential from 
the upper vortex and its image, and similarly the last term contains the contributions from 
the lower vortex and its image. As can be inferred from Fig. [H a necessary condition for 
a steady configuration to exist is that the upper (lower) vortex be of negative (positive) 
circulation, hence only the case T < is of interest to us here. 
In dimensionless variables 
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where the prime notation has been dropped. According to standard theory of point vortices 
in an inviscid fluid, any given vortex moves with the velocity of the flow computed at the 



position of that vortex, excluding its own contribution to the flow. It then follows from 
Eq. (J3J) that the velocity ui = (ui,Vi) of the vortex located at z\ is given by 

ui - ivx = 1 2 ~~ iK \ — " i — + i — ) ' ( 4 ) 
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%l ~ y\ ( V\-V2 V\ V\r\ - y 2 



ui = 1 j K 



r\ + r\- 2(xix 2 + ym) r 2 - 1 1 + r\r\ - 2{x x x 2 + j/i2/ 2 ) 



(5a) 



V\ = -2 — j- + k 



r 



r\ + rl - 2{xix 2 + yiy 2 ) r\ - 1 1 + r^rf - 2(xix 2 + 2/1^/2)/ 

(5b) 

where rf = + 2/ 4 2 , z = 1, 2. The velocity u 2 = (u 2 , v 2 ) of the second vortex is obtained by 
simply interchanging the indexes 1 -H- 2 in Eq. (JSJ) and letting n —k. 

III. DYNAMICS ON THE SYMMETRIC SUBSPACE 

It is not difficult to see from Eq. (JSJ) that if the vortices are initially placed at positions 
symmetrically located with respect to the centerline, i.e., z 2 (0) = ^i(O), then this symmetry 
is preserved for all later times, i.e., z 2 {t) = ~z~i(t) for t > 0. In this section, we study the 
dynamics within this invariant symmetric subspace, where the motion of the lower vortex is 
simply the mirror image of that of the upper vortex with respect to the centerline. Symmetry 
can be enforced experimentally by placing a splitter plate behind the cylinder in the center 
plane of the wake^^. 

With x 2 = X\ and y 2 = —y\, Eq. ([S]) reduces to 

1 1 1 

(6a) 
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Here, the subscripts have been dropped with the understanding that in the remainder of the 

section we restrict our attention to the upper vortex. 
A. Hamiltonian dynamics and phase portrait 



As is well known, the equations of motion for point vortices in a two-dimensional inviscid 
flow, first derived by Kirchhoff, can be formulated as a Hamiltonian system^d The dynamics 
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of point vortices in the presence of closed, rigid boundaries was shown by Lin^P. to be also 
Hamiltonian with the same canonical sympletic structure as in the absence of boundaries. 
For a vortex pair placed in a uniform stream around a circular cylinder, the phase space 
is four- dimensional and has a two-dimensional (2D) invariant subspace corresponding to 
symmetric orbits. The Hamiltonian restricted to the 2D symmetric subspace is given by^i 

H(x,y) = y(l-\) ~-ln - y(r2 ~ 1} . (7) 
1 U> U \ r 2 J 2 v /( r 2_ 1 ) 2 + 4y2 ^ 

The corresponding dynamical equations 

. _ dH_ . _ 
X dy ' ^ dx ' 
where dot denotes time derivative, yield Eq. (Q upon identifying (u,v) with (x,y). 

A phase portrait of this Hamiltonian system for k = 45/32 is presented in Fig. [21 where 
the curves shown are (unevenly spaced) level sets of the Hamiltonian (jJJ). [For convenience, 
these curves were obtained from a direct numerical integration of Eq. A detailed 

description of the main features of this phase portrait will be given below, starting with an 
analysis of the various equilibrium points and their stability. The related problem of the 
symmetric "moving Foppl system," where the cylinder advances through the fluid followed 
by the vortex pair, was recently considered by Shashikanth et al.—, but there the phase 
portrait^ is quite different from the one shown in Fig. [2j because of the additional degrees 
of freedom related to the velocity of the moving cylinder. 

B. Equilibrium points 

The equilibrium positions for the vortex are obtained by setting u = v = in Eq. OH]). 
Three types of equilibrium points can be identified. 

1. Foppl equilibria 

The locus of possible equilibrium positions (xo,2/o) for the upper vortex found by Foppl^ 
is the curve 

r 2 Q - 1 = 2r y , (9) 

with corresponding strength 

^(rg + DW-Q* (1Q) 
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FIG. 2. Phase portrait for the symmetric Foppl system with k = 45/32. The isolated black dots 
are the Foppl equilibria. The dashed curves are the stable and unstable branches of the separatrix 
associated with the equilibrium point on the normal line, and the thick solid lines are the homoclinic 
loops of the equilibrium point at infinity; see text. 



Along the Foppl curve (Q, the vortex strength increases with distance from the center of 
the cylinder and diverges linearly for r -> oo. For the equilibrium point on the edge of the 
cylinder (tq — > 1), the strength vanishes. Notice that Eq. ([9]) yields two branches of solution: 
one in which the vortex pair is behind the cylinder (xq > 0) and the other where the vortex 
pair is in front of the cylinder (xq < 0). The former case models the formation of vortex 
eddies behind a cylinder in a uniform stream and was the primary motivation of Foppl's 
original study^. The latter case has attracted far less attention because it is not usually 
observed in experiments. We note, however, that recirculating eddies are observed in front 
of a circular cylinder near a plane boundary when the gap between the cylinder and the 
plane is sufficiently small^. In this context, the Foppl equlibrium upstream of the cylinder 
may eventually be relevant for flows around a half-cylinder placed on a plane wall (or for 
the closely related situation where a splitter plate is attached to the front of the cylinder), 
although we are unaware of specific experiments in this setting. 
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2. Equilibria on the normal line 



This corresponds to the upper vortex being located on the line bisecting the cylinder 
perpendicularly to the incoming flow^, that is, 

x = 0, y = b, b>l, (11) 

with strength 



2(& 2 - 1)(& 2 

K ~ 



(12) 



b{b A + Ab 2 - 1) ' 

As in the Foppl solution, the strength tends to zero when the edge of the cylinder is reached 
(6—7-1) and diverges linearly with distance from the center of the cylinder. At large 
distances, the vortex strength for this equilibrium is about twice that of a Foppl pair located 
at the same distance from the origin. 



3. Equilibrium at infinity 

Equation (jSJ) also yields equilibrium points at the positions 

x = ±oo, l/oo = -. (13) 

To the best of our knowledge, the existence of this additional equilibrium point at infinity was 
not noted before. Its physical origin, however, can be easily understood, as it corresponds 
to the equilibrium configuration for a vortex pair placed in a uniform stream (without the 
cylinder). At points infinitely far from the cylinder, the flow induced by the image system 
(inside the cylinder) becomes negligible and hence a stationary configuration is possible if 
the vortices with given circulation ±k are placed at the appropriate distance (= k) from 
each other. 



C. Stability analysis 

The linear stability analysis of the equilibria described above is presented next, together 
with a discussion of the nonlinear stability of the Foppl equilibrium. 
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1. Foppl equilibria 



Consider a perturbation of the Foppl equilibrium <Q parameterized as: z = zo + Az, 
where Az = £ + irj, with £ and t] being infinitesimal (real) quantities. Linearization of 
Eq. ([6]) then yields the following dynamical system 



(14) 





where the matrix A reads 



x (r 4 -3r 2 + 2) 
A u = -A 22 = -g , (15) 
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Its eigenvalues A are given by 



_ 4r 8 + 5rg + 2r 4 - 5r 2 + 2 

~ 2? Q ' ^ ' 

2x1(4 + r 2 + 2) 

A 21 - „7^2 i i \ • i 1 ') 
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A 2 = I < Qj (lg) 
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for r > 1. The eigenvalues are thus purely imaginary, and not a complex pair with negative 
real part as found by Foppl^. In other words, the Foppl equilibrium is a center and not 
a stable focus. Our equation (I18p agrees with the expression for the eigenvalues of the 
symmetric modes obtained in Ref. [7] from the linearization of the full 4D dynamical system. 
As can be seen from Fig. [2, the Foppl solution is in fact a nonlinearly stable center, meaning 
that when the vortex is displaced from its equilibrium position by a small (but finite) amount, 
it executes a periodic motion around that point, corresponding to the closed orbits in the 
figure. This periodic motion around the Foppl equilibrium has been observed in numerical 
simulations of the model carried out by de Laat and Coene^. Note that since the eigenvalues 
given in Eq. (118j) do not depend explicitly on the coordinate xq, it follows that the two Foppl 
equilibria, downstream and upstream of the cylinder, have identical stability properties, as 
is evident from Fig. [2j This means, in particular, that if vorticity can be generated upstream 
of the cylinder then stationary recirculating eddies could form in front of the cylinder — a 
situation observed, for instance, in flows around a cylinder placed above a plane wall^ 2 .. 



2. Equilibria on the normal line 



Linearization of Eq. ([6]) around the equilibrium point z = ib yields for the matrix A: 

A u = A 22 = 0, (19) 

6 8 + 106 6 - 8b 4 + Ub 2 - 1 
Al2 ~ &3 (& 2_ 1)(6 4 + 46 2_ 1) > ( 2 °) 

_ 2(^-l)(3fe 2 -l) 

" 63 (6 4 + 46 2_ 1) • (21) 
The eigenvalues A of this matrix are determined by 

_ 2 (36 2 - 1) (6 s + 10& 6 - 86 4 + 146 2 - 1) 
A " 66(6* + 46= -1)2 > °' (22) 

which yields a pair of real eigenvalues, A± = ±\/A^. The equilibrium point on the normal 
line is therefore a saddle, having a stable and unstable direction, as is also evident from the 
phase portrait shown in Fig. [2j The eigenvectors w± associated with the eigenvalues A±, 
respectively, read 



w ± = I ± ^f^ ] . (23) 

Although it was known from numerical simulations^ that the equilibrium point on the nor- 
mal line is unstable (against generic symmetric perturbations), it seems that an explicit 
linear stability analysis for this case was not carried out before, perhaps because these equi- 
libria were not considered physically relevant since they are not observed in experiment^. 
However, when the full nonlinear dynamics is considered, the stable and unstable eigendi- 
rections w± give origin to the respective stable and unstable separatrices, indicated by the 
dashed curves in Fig. [2J In this sense, the existence of an equilibrium point on the normal 
line is dynamically felt by a vortex even if it is placed far from this "unphysical" equilibrium. 

3. Equilibrium at infinity 

The matrix A of the linearized system around the equilibrium point at infinity is given 

by 

9/011 

(24) 
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FIG. 3. Trajectories near a nilpotent saddle. 



which is nilpotent and has two zero eigenvalues. To study the stability of this equilibrium 
point, one needs to examine the nonlinear contributions. To this end, we note that for 
|x| — > oo and y ~ y^, Eq. (16bl) assumes the form 

* = -?■ (25) 

It then follows from a theorem in ordinary differential equations^ that, in view of the cubic 
term in Eq. ( 1251) . the equilibrium point is a degenerate or nilpotent saddle^, for which the 
two eigenvectors are the same. The behavior of trajectories in the neighborhood of a generic 
nilpotent saddle is illustrated in Fig. [3j The behavior near the nilpotent saddle at x = ±00 
and y = y^ can be described as follows. A vortex placed very far downstream and below 
(above) the line y = y^ will move away from (towards) the equilibrium point at x = 00. 
Similarly, a vortex placed very far upstream will move away from (towards) the equilibrium 
point at x = -00 if y > y^ (y < y^). 

The stable and unstable separatrices associated with the nilpotent saddle at infinity form 
two homoclinic loops^, called nilpotent saddle loops, which are indicated in Fig. [2] by thick 
solid lines and correspond to the level curves passing through this equilibrium point: 



The nilpotent saddle loops encircle the Foppl equilibria and define their region of nonlinear 
stability, in the sense that vortex trajectories are closed for initial positions inside the loops 




(26) 
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and unbounded otherwise. In this way, the nilpotent saddle at infinity, which went unnoticed 
until now, allows us to fully characterize the nonlinear stability of the Foppl equilibrium. 

For unbounded orbits, the long-time asymptotic behavior depends on the location of the 
vortex initial position with respect to the separatrices associated with the equilibrium point 
on the normal line. A vortex placed downstream of the cylinder between the nilpotent 
saddle loop and the separatrices of the equilibrium point on normal line will eventually be 
convected away by the free stream; see Fig. [2j In particular, if the vortex starts very far 
behind the cylinder at a position that is below the nilpotent saddle loop and above the 
stable separatrix, it first moves towards the cylinder, turns around the Foppl equilibrium, 
and is then "reflected" back to infinity. Even more surprising trajectories arise if the vortex 
is placed downstream below the stable separatrix, for it will be close enough to its image 
below the centerline to be able to overcome the cylinder and move off to infinity upstream. 
(A related phenomenon occurs in the inviscid coupled motion of a cylinder initially at rest 
and a vortex pair starting at infinity with no imposed background flow^. When the cylinder 
is less dense than the fluid, it is found that if the vortices are released sufficiently above the 
centerline they reverse relative to the moving cylinder; otherwise, they move over and past 
the cylinder.) Unbounded trajectories for the Foppl system also result for initial positions 
upstream of the cylinder: i) if placed above the stable separatrix, the vortex moves down- 
stream to infinity; and ii) if placed between the stable separatrix and the nilpotent saddle 
loop, the vortex goes around the Foppl equlibrium in front of the cylinder and returns to 
infinity upstream; see Fig. |2j It is again the hitherto unnoticed nilpotent saddle at infinity, 
together with the precise nature of the equilibrium point on the normal line, that allows us 
to go beyond linear stability analysis and capture the full phase portrait in the symmetric 
subspace. 

We stress that closed orbits exist only when the flow is symmetric. Nonsymmetric per- 
turbations inevitably cause the vortex pair to move off to infinity, as we demonstrate next. 

IV. NONSYMMETRIC DYNAMICS 

In this section, the effect of antisymmetric perturbations on the equilibria of the Foppl 
system is studied. We begin by observing that the dynamics of two counter-rotating point 
vortices possesses a conjugation symmetry. To describe this symmetry, let Zi(t] z lt0 , z 2 ,o) an d 
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Z2{t] £1,0 > z 2,o) denote the trajectories of the upper and lower vortices, respectively, with initial 
positions Z\$ and £2,0- For the dynamical system defined by Eq. (jSJ) and the corresponding 
equation for the second vortex, one can verify that the following relations hold 

Zi{t; ^2,0) £1,0) = z%(t] zifi, Z2,o), (27a) 

Z2{t; Z2,o, zifl) = zi(t) zifi, Z2,o). (27b) 

In other words, for any given pair of initial positions, z lt0 and z 2 ^, there exists a "conjugate 
pair" of initial positions, 2:2,0 and ~Zi,o, such that the vortex trajectories of the first pair are 
the complex conjugate of those of the second pair. 

Any perturbation of a vortex-pair equilibrium can be written as the superposition of a 
symmetric perturbation and an antisymmetric one. To be precise, antisymmetric perturba- 
tions are of the form 

z\ = z + Az, z 2 = z Q - Az, (28) 

where z denotes a generic equilibrium point and Az — £ + irj. Since the antisymmetric 
subspace of the full 4D phase space is invariant under linear dynamics, we can focus on the 
upper vortex in carrying out our linear stability analysis. 



A. Foppl equilibria 

Linearization of Eq. (JSJ) around the Foppl equilibrium ()9]) with respect to antisymmetric 
perturbations (j2"8~j) yields 



V J \ V 



where the matrix B is given by 



Bn = -B 22 = ,30) 
_ 3r 6 - 5rg + 2 

~ ^9 ' ( 31 ) 

4r 8 + 3r 6 - 4r 4 - 5r 2 + 2 

B 21 = — . (32) 





This matrix has a pair of real eigenvalues, A± = ±y/X*, where 



A 2 = *£±Mz±i±l. (33) 
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FIG. 4. Vortex trajectories for antisymmetric perturbations of the Foppl equilibrium for k = 45/32, 
in which case xq = y/55/4 and yo = ±3/4 (black dots). The solid and dashed curves are the 
trajectories starting along the unstable directions w + and — w+, respectively, while the short 
straight lines indicate the axes defined by the stable direction w_. The dotted lines represent the 
loci of the Foppl equilibria. 

The Foppl equilibrium is therefore a saddle with respect to antisymmetric perturbations, 
while it is a center with respect to symmetric perturbations, as seen earlier. That is, the 
Foppl equilibrium is a saddle-center of the full 4D dynamical system^. We note in passing 
that, although Foppl obtained a pair of real eigenvalues for the case of antisymmetric per- 
turbations, his original formulae for the eigenvalues are in error—. Our expression (1331) is 
in agreement with the eigenvalues of the skew-symmetric modes obtained by Smittf from 
the linearization of the full dynamical system. The eigenvectors w± associated with the 
eigenvalues A± are readily computed, with the result 

w± = I ^ . (34) 

In Fig. HI we show in solid curves the pair of vortex trajectories obtained by slightly 
displacing the vortices from their equilibrium positions in the directions defined by the 
unstable eigenvector w + , while the trajectories obtained by slightly displacing the vortices 
in the opposite directions are shown in dashed curves. The latter pair of trajectories is the 
complex conjugate of the former by conjugation symmetry. Note that for the first pair of 
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trajectories, the lower vortex initially moves towards the centerline and upstream, while the 
upper vortex moves away from the centerline and downstream. At later times, the vortex 
pair moves off to infinity with the lower vortex trailing behind the upper vortex. For the 
second pair of trajectories, the upper and lower vortices switch role; see Fig. |H In the flow 
of a real fluid past a cylinder, the two basic instabilities associated with displacements along 
the unstable directions ±w + happen alternately and constitute the origin of vortex shedding 
that leads to the formation of the Karman vortex street^. In like manner, the suppression of 
vortex shedding by placing a splitter plate behind the cylinder— ^ is consistent with the fact 
that the Foppl equlibrium is nonlinearly stable with respect to symmetric perturbations; see 
Sec. [V] for further discussions on vortex shedding and its suppression by a splitter plate. 

For small, generic antisymmetric perturbations, the vortices move along trajectories that 
follow closely the ones depicted in Fig. HI Whether a vortex pair eventually moves up or down 
is determined by the initial position of the upper vortex relative to the stable direction w_, 
which is indicated in Fig. H] by the short straight line passing through the Foppl equilibrium. 
If the initial position of the upper vortex is to the right (left) of the stable direction, then the 
vortex pair asymptotically moves upwards (downwards). This explains the behavior seen in 
the numerical simulations reported in Ref. [25], where nearby initial positions around the 
Foppl equilibrium were found to lead to close-by trajectories. 

Since any degree of antisymmetry in the initial perturbation causes the vortex pair to 
move off to infinity, the Foppl equilibrium is unstable under generic perturbations. As an 
example, Fig. [5] shows vortex trajectories obtained by displacing the Foppl pair (at r = 2) 
by the amounts Az\ = A22 = — 0.25 + iO. 005. During the linear stage, the trajectories are a 
superposition of a symmetric orbit and a growing mode associated with the antisymmetric 
component of the perturbation, which ultimately leads to asymptotic trajectories with the 
vortices moving parallel to each other. 

B. Equilibria on the normal line 

For antisymmetric perturbations of the equilibrium ( 1TTT) on the normal line, the matrix 
B assumes the form 

B u = B 22 = 0, (35) 
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FIG. 5. Trajectories resulting from a generic perturbation Azi = A^2 = —0.25 + i0.005 of the 
Foppl pair at ro = 2 (black dots). 

2(3& 6 + 6 4 + 56 2_ 1) 
12 63( fe 2_ 1 )( 6 4 + 46 2_ 1 )' 

6 2 - 1 

5 21 = — (37) 

with eigenvalues A given by 

_ 2 (3fe 6 + + 56 2 - 1) 
This yields a pair of real eigenvalues, A± = i-s/A 2 , with respective eigenvectors: 

(39) 

In Fig. [61 we show the vortex trajectories (solid curves) obtained by slightly displacing 
the vortices from their equilibrium position along the unstable direction w + for 6 = 2. The 
initial motion here is somewhat similar to what is seen for a Foppl pair, in the sense that 
one vortex moves upstream towards the centerline and the other moves downstream away 
from the centerline. The main difference is that for later times, the vortices now end up 
moving upstream. The long-time dynamics in this case is also more sensitive on the initial 
conditions: for somewhat larger perturbations, the vortices are eventually carried away by 
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FIG. 6. Vortex trajectories (solid curves) associated with the unstable direction w + of the equilib- 
rium at z = ±2i (black dots). The dashed curves are trajectories resulting from the antisymmetric 
perturbation Az = 0.16. 

the free stream. An example where this happens is indicated by the dashed curves in Fig. El 
which represent the vortex trajectories for the antisymmetric perturbation Az = 0.16. 

As already argued in Sec. IIII C 2\ although the equilibrium point on the normal line is 
not directly observed in experiments, it is important to know its instability properties under 
both symmetric and antisymmetric perturbations. This knowledge contributes to a better 
understanding not only of the full nonlinear dynamics of the Foppl system but also of more 
general flows, such as the case of stationary vortex patches above and below the cylinder in 
a uniform stream, where similar unstable modes are observed^. 

V. DISCUSSION AND CONCLUSIONS 

In this paper, we have investigated a two-dimensional vortex model for the formation of 
recirculating eddies behind a fixed cylinder placed on a uniform stream. The model, which 
was first introduced by Foppl^ almost a century ago, has two main simplifying assumptions: 
i) the fluid is treated as inviscid and hence the flow is potential, and ii) the size of the 
vortex core is neglected and so the vortices are considered to be point-like. In spite of 
these simplifications, the model is known to be in qualitative agreement with real flows 

17 




past a cylinder, as was already pointed out by Foppl in his original paper. Several novel 
features of the Foppl model have been obtained in the present work, which help one to 
better understand the basic dynamics of vortex shedding behind a cylinder. 

In real flows, governed by the Navier-Stokes equations, stationary vortices behind a cylin- 
der are formed at moderate Reynolds number (Re < 50). As the Reynolds number increases 
past Re ~ 50, the configuration loses its symmetry and becomes unstable. New vortices 
then start to form alternately on both sides of the cylinder, while the vortices further down- 
stream break away and develop into a Karman vortex street, as described by Foppl^. It has 
been argued by Roshkc3 that "possibly the breaking away should be regarded as primary, 
resulting in asymmetry." The analysis presented in Sec. IIVA1 makes it clear that the reverse 
scenario is more plausible: the asymmetrical disturbances induce the instability of the vortex 
pair which then breaks away from the cylinder. As vorticity is continuously generated from 
the separated boundary layer on both sides of the cylinder, new vortices are formed and 
alternately shed into the far wake of the cylinder according to the unstable modes shown 
in Fig. |H Direct numerical simulations (DNS) of two-dimensional flows past a cylinder per- 
formed by Tang and Aubry^ have confirmed that the mechanism for the instability of the 
symmetric eddies in real flows is qualitatively described by the instability of the point- vortex 
model. 

It is experimentally observecr^^ that vortex shedding is suppressed if a splitter plate is 
installed behind the cylinder in the center plane of the wake. The presence of the splitter 
plate tends to enforce symmetry of the flow with respect to the centerline, thus effectively 
reducing the appearance of antisymmetric disturbances behind the cylinder. The suppres- 
sion of vortex shedding in this case is thus entirely consistent with the fact that the Foppl 
equilibria of the vortex-point model is nonlinearly stable against symmetric perturbations 
and that vortex shedding is induced by unstable antisymmetric modes, as discussed above. 
This scenario has been confirmed by DNS of flows past a cylinder with symmetry imposed 
along the centerline recently performed by Kumar et al.—. The problem of stationary con- 
figurations for vortex flows past a cylinder with patches of constant vorticity has also been 
studied numerically by Elcrat et al.—^Z-. These authors found two families of solutions, rep- 
resenting desingularized versions of the Foppl and the normal equilibria, respectively, which 
have the same stability properties as the corresponding point-vortex equilibria. 

In conclusion, we have seen that the Foppl model, where a pair of counter-rotating point 
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vortices move around a circular cylinder in the presence of a uniform stream, is a rich 
nonlinear dynamical system whose features — notably its stability properties — bear a direct 
relevance to our understanding of the vortex shedding mechanism in real flows. The results 
obtained here should, in principle, carry over to more general geometries, such as vortex 
motion around a plate or around a cylinder with noncircular cross section. 
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